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Abstract 

A two-dimensional earthquake model that consists of a single block resting upon a 
slowly moving rough surface and connected by two springs to rigid supports is stud- 
ied. Depending on the elastic anisotropy and the friction force three generic regimes 
are possible: i) pure creep; ii) pure stick-slip motion; and iii) a mixed regime. In all 
cases the long-time dynamics (fixed point, periodic orbit or chaos) is determined by 
the direction of the pulling velocity. The possible relevance of our findings to real 
faults is briefly discussed. 
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1 Introduction 



Since the seminal work of Burridge and Knopoff [1], spring-block models have 
been recognized as useful tools to study earthquake dynamics [2] . Because they 
are simple to treat, both theoretically and computationally, spring-block mod- 
els can provide important physical insights that might otherwise be much more 
difficult to obtain. Models with many blocks have been used to investigate 
the origin of certain power laws that appear in the statistics of earthquakes, 
such as the Gutenberg- Righter law for the size-distribution of earthquakes [2]. 
Spring-block models with only a few degrees of freedom, on the other hand, 
can probe the basic dynamics of frictional sliding and are also of interest in 
their own right as examples of low- dimensional nonlinear dynamical systems. 
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For instance, two-block models have been found to exhibit chaos [3,4,5], thus 
suggesting that actual faulting might perhaps be a chaotic phenomena. Even 
in the simplest case of a one-block model (with the block moving on a line) 
the dynamics is highly nontrivial. Indeed, this system displays a discontinuous 
transition from creep to stick-slip motion as one varies a parameter governing 
the dependency of the friction force on the velocity [6]. A two-dimensional 
version of the one-block model was recently considered [7] and a richer dy- 
namics was found, including a chaotic regime where stick-slip events occur 
intermittently amidst creep. 

In this paper we present a detailed study of a 2D one-block model for earth- 
quakes, in which the block rests upon a moving planar surface and is connected 
by springs to two perpendicular rigid walls, with a velocity-weakening friction 
force acting between the slider and the surface. An earlier discussion on this 
model was reported by Ryabov and Ito [7]. Here we give a much more com- 
plete analysis and present several novel and surprising results. We show that 
in the limit of very slow pulling the model is governed by three dimensionless 
parameters, namely, the asymmetry parameter k corresponding to the ratio 
between the two spring constants, a parameter 7 governing the decrease of 
friction with velocity, and the direction 6y of the pulling velocity. Depending 
on the values of k and 7, the model displays three possible generic regimes: 
i) pure creep; ii) pure stick-slip motion; and ill) a mixed regime where both 
creep and stick-slip occur. In each of these cases, the long-time dynamics (fixed 
point, periodic orbit or chaos) is starkly dependent on Ov, as we will see below. 



2 The model 



The model we consider consists of a block of mass m connected to motionless 
walls by two springs of stiffness and ky^ as shown in Fig. la. The position 
of the block is described in a system of coordinates (x, y) fixed to the spring 
supports, with the origin being placed at the point where the elastic force 
vanishes. The block rests upon a surface that moves with a constant velocity 

— * — * _ 

V, where \V\ <^ 1. There is friction between the block and the moving surface, 
so that initially the block moves with the substrate until the elastic force F^i 
overcomes the static friction force Fq, at which point the block slips with 
respect to the surface. In the slip phase of the motion, we adopt a commonly 
used velocity- weakening friction law [8] , in which the magnitude of the friction 

force Ffr is given by |Ffr| = Fo^{vr/Vi), where — \Jii — V^)^ + iv — VyY is 
the relative velocity of the slider with respect to the surface, Vf is a typical 
velocity scale for the friction force, and ^{x) = + x). 

In the remainder of the paper we work with non-dimensional quantities, where 
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Fig. 1. (a) Two-dimensional spring-block model and (b) the ellipse 
C = {{x,y) I K^x'^ + = 1} where the elastic force equals the static fric- 
tion. Starting in a point inside C, the block is dragged by the surface along the 
direction O^, of the pulling velocity until it hits a point on C, whereupon slip begins. 
The angle 9 denotes the direction opposite to the elastic force at this point. 



we have rescaled the length and time scales by Fo/ky and Jm/ky, respectively, 



and velocities by Fo/ymky. In these units, the dimensionless parameters char- 

acterizing the problem are the pulling velocity z7, the anisotropy parameter 
K, = kx/ky, and the friction parameter 7 = l/2i/f, where z/f is the dimension- 
less characteristic velocity of the friction. If we now denote by C the curve of 
maximum displacement, where the condition |Fei| = Fq is verified, it then fol- 
lows that C corresponds (in dimensionless units) to the ellipse k^x"^ + = 1; 
see Fig. lb. Starting from relative rest at a point inside C, the block will be 
dragged by the surface until it reaches a point on C, upon which the block 
starts to slide. In the slip phase the equations of motion are given by 

X = — + $(27fr) C0S(y9fr, (1) 




where the angle gives the direction of the friction force. These equations 
represent two coupled (nonlinear) oscillators, and in this sense the system 
is somewhat similar to the two-block model considered in Ref. [3]. Our 2D 
one-block model is, however, qualitatively richer in that it contains a new 
parameter, namely, the direction of the pulling velocity which strongly affects 
the dynamics of the system, as we will see shortly. 

A final caveat about the model concerns the choice of direction of the friction 
force. The problem comes about because at the onset of motion the friction 
force is opposite to the elastic force, whereas during the slip phase of the 
motion the friction force is expected to point in the direction opposite to the 
relative velocity. Since the elastic force is not in general collinear with the 
velocity, we should thus prescribe a mechanism that allows the switching of 
direction of the friction force. Here we follow Ryabov and Ito [7] and assume 







y^-y + $(27Wr)sin9?fi 



(2) 
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that in the acceleration phase of the block motion, the direction ipf^ of the 
friction force is chosen according to the following rule 



</^/r- = TT + ^Pel - {(Pel " Vvr)[i " exp{-QVyr)] 



(3) 



where (p^r ^-nd ipei denote the directions of the relative velocity and the elastic 
force, respectively, and is a parameter describing the 'memory' of the friction 
force. (The qualitative behavior of the model does not depend on the value 
of fl and so we have used i7 = 50 throughout the paper.) In the deceleration 
phase, the friction force is kept opposite to the velocity to ensure that the 
block always comes to a stop. 

The dynamics of the model described above can be conveniently described in 
terms of a one- dimensional map. To see this, let us first parametrize the curve 
C in terms of the angle 9 defined by the direction opposite to the elastic force 
with respect to the x axis; see Fig. lb. As a function of 9, the coordinates 
(xq, Uq) of a point on C are given by the following parametric equations: Xq — 
K,~^cos9, Uq = sin9. Now suppose that the block starts to slip from a given 
point on C labeled by 9. The block will slide for a while until it eventually 
comes to a stop with respect to the moving surface. After this, the block 
sticks to the surface and is brought back to another point 9' on C, where a 
new slip cycle begins, and so on. Thus, the overall dynamics of the model can 
be described by a one-dimensional map 9' = f{9). Obviously the mapping 
function f{9) cannot be computed explicitly since this entails solving (1) and 
(2) during the slip phase of the motion. However, for a given set of parameters 
K, 7, and P, the map f{9) can be easily constructed on the computer [7]. 



3 Analysis of the model and results 

We begin the analysis of our model by considering the early stages of a slip 
event. Suppose that the block, as it is being dragged by the moving surface, 
reaches a point (xo,l/o) G C where relative motion begins. At the onset of slip 
the block will move in the direction of the elastic force F^i = {—KXo,—yo), 
while the friction force will act in the opposite direction. In order to study 
the nature of the ensuing motion, we linearize the equations of motion in the 
directions parallel and perpendicular to the elastic force. To accomplish this, 
it is convenient first to switch to the reference frame where the substrate is at 
rest and then write the equations of motion in the coordinates {i,r]), defined 
as the block displacement in the directions parallel and perpendicular to the 
elastic force, respectively. The resulting linearized equations are 



(4) 
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Fig. 2. The ellipse C of maximum displacement (solid line) and the critical ellipse 
(dashed line). We have 7^ < 1 in (a) and 1 < 7^ < k in (b). 

fj = K{n- l)xoyoi - f^-if^-xl + yl)v + Pt: (5) 
where 



A = K^xl + yQ, a = -{n'^xou^ + yoiyy), = K{yoVx - yo^y)- (6) 



We see from (5) that in the limit i/ — > the transversal coordinate 77 can grow 
only if so does the component ^. We can thus set 77 = in (4), so that the 
linear equation for the variable ^ becomes 

-(-2^i + Ai = at. (7) 



Equation (7) for the case A = 1 was studied in great detail by one of the 
present authors [6]. There, it was shown that in the limit that i/ — > the 
system undergoes a phase transition at 7 = 1 in the following sense: for 7 < 1 
there is only creep, meaning that the linear approximation is always valid and 
the amplitude of a slip event vanishes as z/ — 0, whereas for 7 > 1 one has 
stick-slip motion in which case the slip amplitude remains finite as — > 0. 
These results can be trivially extended to the general case given in (7). Here 
the critical point is reached when the condition 7^ = ^4 is satisfied, so that 
if 7^ < A we have creep, while for 7^ > A stick-slip occurs. Note, however, 
that since the parameter A depends on the point (a;o,yo) where the slip was 
initiated, it follows that for a critical point to exist the ellipse defined by the 
critical condition A = 7^ must intersect the curve C. If such an intersection 
occurs, then the model displays both creep and stick-slip, otherwise the motion 
is either pure creep or pure stick-slip, as described next. 

Pure creep: 7 < 1. In this case the critical ellipse /I = 7^ is inside C (Fig. 2a), 
so that ^"^ < A for any point (xq, yo) G C, and hence the linear approximation 
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Fig. 3. Bifurcation diagram for the map f{9) in terms of the control parameter 9i,. 
In (a) we have 7 = 6.5, k = 50, u = 10"'', corresponding to the creep and stick-slip 
regime, and the dashed line denotes the critical point 6c = 0.41. The inset shows 
the initial part of the bifurcation diagram, with the solid straight line representing 
the 45° line. In (b) we have stick-slip only with 7 = 6.5, k = 18.95, and u = 10~^. 

(7) always holds (in the limit — > 0). The long-time dynamics of the system 
is quite simple and can be easily described in terms of the map 6' = f{6) 
introduced above. Let us define by 9„ the direction of the pulling velocity. Then 
for any initial condition the block will reach the fixed point 9 — 9* — O^,, where 
the elastic force is opposite to V (see Fig. 2a), so that the resulting motion is 
essentially one-dimensional creep: starting from the point on C labeled by 9*, 
the block slides a very small distance along the direction of the elastic force, 
stops, and then is brought back to the same point 9* , and so on. 

Creep and stick-slip: 1 < 7^ < k. In this situation the critical ellipse and 
C intersect each other at four critical points 91, where tan^ = ^5^; see 
Fig. 2b. Without loss of generality, let us assume that 9y is restricted to the 
first quadrant, so that for < 6* < 6*1 we have creep (7^ < A), whereas for 
9{ < 9 < 7r/2 stick-slip occurs (7^ > A). Thus, for 9^ e [0,9c) the creep 
fixed-point 9 = 9* = 9i, is the only attractor. This is illustrated in the inset 
of Fig. 3a, where we plot the initial part of the bifurcation diagram for the 
map f{9) with 7 = 6.5, k = 50, u = 10^^, and 9^ as the control parameter. 
In this inset, the dashed line represents the critical point 9c = 0.41 and we 
see that for 9i, < 9c the fixed point indeed follows closely the 45° hne (solid 
straight line), with the small deviation being caused by the finite value of u 
that smoothes out the transition from creep to stick-slip. As 9iy increases past 
9c, the system remains in a stick-slip fixed point for awhile but eventually 
undergoes a sequence of period-doubling bifurcations leading to chaos, as seen 
in the main bifurcation diagram shown in Fig. 3a. For larger 9,^, the system 
behaves intermittently [7] in the following sense: the block spends a long time 
in the creep region (below the dashed line) until it eventually crosses the 
critical point 9c, after which it undergoes a large slip event and is reinjected 
back into the creep region, and so on. Finally, as 9i, approaches 71/ 2 there is a 
reverse period-doubling cascade back to a fixed point. 
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Fig. 4. Phase diagram in the plane {6^, k) for the pure stick-slip regime. FP denotes 
the region with a stable fixed point, Pn indicates regions of period-n orbits, and 
the shaded areas represent chaotic regions. Here we used 7 = 6.5 and u = 10"^. 



Pure stick-slip: 7^ > k. Here the critical ellipse is exterior to C and so we 
have 7^ > A for any point on C, hence the block always undergoes stick-slip 
motion. As the parameters k, and 61, are varied (with 7 fixed), the system now 
displays a very rich behavior that is summarized in the phase diagram shown in 
Fig. 4, whose main features we shall briefly discuss. First we see that for small 
values of 9^, the system always reaches a stable fixed point (region labeled by 
FP), where slip events occur periodically. As 9i, increases (with k fixed), what 
follows next depends on the value of k. For k very close to 1 the fixed point 
remains stable for all values of O,^; sec bottom of Fig. 4. However, for larger k 
the fixed point eventually goes unstable and a orbit of period two is born (re- 
gion P2). The fate of this period-2 orbit, as 9^, increases further, also depends 
on the value of k. For small k, e.g., k = 10, the system simply undergoes a 
reverse period-doubling bifurcation back to the fixed point. For intermediate 
values of 16, one observes a full period-doubling cascade leading to 

chaos (shaded region), followed by a reverse cascade. For larger k, windows of 
periodic orbits and their associated cascades appear inside the chaotic region. 
Several of these windows are clearly visible in Fig. 4, although we have labeled 
only the first three of them, namely, the windows of period 3, 4, and 5. Note 
also that depending on the value of k the secondary period-doubling cascades 
may not develop fully. An example of this is given in Fig. 3b where we show 
the bifurcation diagram for the map f{9) with 7 — 6.5 and k — 18.95. 
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4 Conclusions 



In conclusion, we have studied a two-dimensional one-block model for earth- 
quakes that shows a very rich behavior and is able to reproduce a host of 
relevant dynamical regimes such as creep, stick-slip (periodic and chaotic) 
and intermittent stick-slip amidst creep. Particularly interesting is our finding 
that the dynamics of the model is strongly dependent on the direction of the 
pulling velocity. This property thus suggests that, in addition to the friction 
force acting on a fault and the stiffness of the loading system, the direction 
of shear might also play an important role in determining the fault seismic 
activity (or lack thereof). Such mechanism could perhaps help explaining, for 
instance, why different segments of a fault system may exhibit distinct seismic 
patterns although they presumably have similar geophysical properties. These 
are interesting possibilities that certainly deserve to be investigated further. 
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